Multifractal burst in the spatio-temporal dynamics of jerky flow 
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The collective behavior of dislocations in jerky flow is studied in Al-Mg polycrystalline samples 
subjected to constant strain rate tests. Complementary dynamical, statistical and multifractal 
analyses are carried out on the stress-time series recorded during jerky flow to characterize the 
distinct spatio-temporal dynamical regimes. It is shown that the hopping type B and the propagating 
type A bands correspond to chaotic and self-organized critical states respectively. The crossover 
between these types of bands is identified by a large spread in the multifractal spectrum. These 
results are interpreted on the basis of competing scales and mechanisms. 

PACS numbers: 62.20Fe, 05.45Df, 05.45Tp, 0.54Ac 



Jerky flow or the repeated yielding of alloys during 
plastic deformation has been studied since the turn of 
the last century and still continues to engage the at- 
tention of scientists. The Portevin - Le Chatelier effect 
(PLC), as it is referred to, has been observed in many di- 
lute alloys [^| . It is one of the few striking examples of the 
complexity of the spatio-temporal dynamics arising from 
the collective behavior of defect populations. In a tension 
test with a constant imposed strain rate e a , in practice 
a constant pulling velocity V, the effect manifests itself 
as a series of serrations in the stress-time or strain curve. 
Each stress drop is associated with the nucleation of a 
band of localized plastic deformation which under certain 
conditions propagates along the sample. The continued 
interest in the PLC effect arises not only from its intrigu- 
ing spatio-temporal behavior but also from its detrimen- 
tal influence on the mechanical properties of materials. 
Our purpose is to quantify the connection between the 
spatial (types of deformation bands) and temporal (na- 
ture of serrations) manifestations of the phenomenon. 

It is now accepted || that the microscopic origin of 
the PLC effect is the dynamic aging of the material due to 
the interaction between mobile dislocations and diffusing 
solute atoms. In a certain range of strain rates and tem- 
peratures, the diffusion time of the solute atoms is of the 
order of the waiting time of the dislocations temporar- 
ily arrested at obstacles. Then, solute atoms can diffuse 
to, and age the dislocations. Thus, the stress needed 
to unpin the dislocations increases with increasing wait- 
ing time and decreasing strain rate. At the macroscopic 
scale, this inverse force versus flux relation translates into 
a negative sensitivity of the flow stress to the strain rate. 
When the imposed strain rate falls into this range, the 
plastic deformation becomes nonuniform. 

In polycrystals, the metallurgical taxonomy distin- 
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guishes three generic types of serrations, namely types A, 
B and C Q . On increasing strain rates or decreasing the 
temperature, one first finds the type C, identified with 
randomly nucleated static bands and large characteristic 
stress drops. Then the type B with smaller serration am- 
plitudes is found. The bands formed are still localized 
and static in nature, but forming ahead of the previous 
band in a spatially correlated way giving the visual im- 
pression of a hopping propagation. Finally, one observes 
continuously propagating type A bands associated with 
small stress drops. The decrease in the amplitude of the 
stress drop with increasing e a reflects the decrease in the 
solute concentration on the arrested dislocations. 

This wealth of spatio-temporal features has long defied 
a proper understanding. Renewed interest has come from 
the connections made to nonlinear dynamics. Dynamical 
|| ||, fjj and statistical methods ||, fj| 0, |9| have recently 
been applied to characterize the complexity of the stress - 
time series. Studies in Refs. |IJ g, ||, fj], ^ suggest that two 
distinct dynamical regimes encountered in single crystals 
might also be found in polycrystals. One is the chaotic 
dynamics |Q characterized by only a few degrees of free- 
dom, identified as stress and dislocation densities [nO 
The other is self-organized critical (SOC) dynamics 
involving infinite degrees of freedom ]Tl[ | . In polycrystals 
only chaotic dynamics has been identified so far ||. 

This paper is focused on type B and A bands in poly- 
crystalline samples as they are expected to show interest- 
ing dynamical behavior in contrast to the uncorrelated 
type C bands not considered here. We propose to iden- 
tify the different dynamical regimes and their correla- 
tions with type B and A bands in polycrystalline sam- 
ples. To unravel this missing connection, dynamical, sta- 
tistical and multifractal analyses of the stress-time series 
are combined with the identification of the spatial band 
patterns. The focus is on the multifractal analysis. The 
motivation for using this method stems from a conceptual 
similarity between the present transition from hopping to 
propagating bands and the Anderson transition in disor- 



2 



dered systems (T^]. In the Anderson model JL3|, wave 
functions are localized when the energy, E, is below the 
mobility edge E c and extended for E > E c . At E = E c , 
the states are shown to exhibit a multifractal character 
|0. In the PLC effect, the hopping type B bands are 
essentially localized, whereas type A bands are delocalized 
in the sense that they are propagating. It is this crossover 
from localized to delocalized nature of the bands that we 
wish to capture by multifractal analysis. We identify the 
dynamical regimes associated with types B and A bands 
as chaotic and SOC type respectively, and show that the 
range of multifractality exhibits a sharp peak in the tran- 
sition region. This is shown to be a consequence of a few 
relevant competing scales and mechanisms. 

The tensile samples of an Al — 2.5%Mg alloy, of length 
L = 70 mm, were cut out of a polycrystalline cold-rolled 
sheet, parallel to the rolling direction. The grains were 
anisotropic in shape, with an aspect ratio of 5. Tests 
were carried out at room temperature for eight values 
of e a — V/L in the range of type B and A bands, from 
5.56 x 10~ 6 s -1 to 1.4 x 10 _2 s _1 . The stress was recorded 
at a sampling rate of 20 Hz - 200 Hz. Three typical stress- 
time, o~(t), curves (corrected for the drift due to strain 
hardening ) are shown in Fig. 1 together with \da/dt\. 
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Stress-time series at three strain rates: 



(a) 5.56xl0 _6 s _1 , (c)2.8 x lO^s" 1 and (e) 5.56 x 10 _3 s _1 . 

(b) ,(d) and (f) are the corresponding plots of \da/dt\. 

The dynamical analysis starts by embedding the time 
series in a higher dimensional space using the time de- 
lay technique. Consider the stress signal {a(k),k = 
1, 2, K}, where K is the number of data points and 
k is in units of the sampling time. Then, a d-dimensional 
vector is defined by = [<r(k),a(k + r), cr(fc + (d — 
l)r)],fc = 1, (K — (d — 1)t), where r is a delay 
time (in units of the sampling time) . The set , k — 
1, (K — (d — 1)t)} constitutes the reconstructed at- 
tractor. The correlation dimension, v, is obtained by 
using the Grassberger-Procaccia algorithm pjfl . In this 
method one calculates the correlation integral C(r), de- 



fined as the fraction of the pairs of vectors (£, , £j ) closer 
than a specified value r. The self similar structure of the 
attractor, when it exists, is revealed by the scaling rela- 
tion C(r) ~ r v in the limit of small r. As the embedding 
dimension d is increased, the slope lnC(r)/lnr tends to 
a constant value taken as the correlation dimension. The 
Lyapunov spectrum is computed using Eckmann's algo- 
rithm [fL5| , suitably modified for short noisy time series. 
In our algorithm, the number of neighbors of a vector £j 
contained in a shell centered on £j (whose size e s is mea- 
sured in % of the size of the attractor) is large enough to 
properly sample the statistics of uncorrelated noise cor- 
rupting the original signal. We consider the dynamics to 
be deterministic if a fair range of e s can be found such 
that stable values (i.e. constant in that range) emerge 
for a positive exponent and a zero exponent M . 



In the region of strain rates 5.6 x 10 s < e a < 1.4 x 
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1 , the data sets typically contain 10000 to 12000 
The band patterns are of type B. We illustrate 
the results with the data file at e a — 5.6 x 10 _6 s _1 . Fig. 
2 shows the log-log plot of C(r) for d = 15 to 18 using 
t = 8. The slopes are seen to converge to a value v ~ 4.6 
for d = 17 and 18 in the range —4.9 < Inr < —3.2. 
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FIG. 2: Correlation integral C(r) using d — 15 to 18, r = 8 for 
e a = 5.56 x 10 _6 s _1 . The curves corresponding to d = 15 to 
17 have been displaced with respect to d = 18 by a constant 
amount. Dashed lines are guide to the eye. 

The Lyapunov spectrum shown in Fig. 3 has been cal- 
culated using d = 5. Stable positive and zero Lyapunov 
exponents are seen in the range 6% < e s < 12%. The 
Lyapunov dimension Dry obtained from the spectrum 
by using the Kaplan- Yorke conjecture turns out to be 
Dky ~ 4.6 = v. Therefore, we conclude that the time 
series is of chaotic origin. Similar results were obtained 
with other samples in this region of e a , with values of 
v and Dky m the range 4.4 to 4.6. In contrast, the 
correlation dimension did not converge and no positive 
Lyapunov exponent could be found for the data sets at 
higher e Q , implying that the dynamics is no more chaotic 
in that region of strain rate. In these cases, only the 
largest Lyapunov exponent could be calculated since the 
data contain only 3000 to 5000 points. 

The quantity \da/dt\ clearly reflects the bursts of plas- 
tic activity (see Fig. 1). Therefore, we use its finite 
difference approximant ipi(U) for statistical and multi- 
fractal analysis. Let Aip denote the amplitude of the 
bursts and At their durations. We investigate the distri- 
butions D(Aip) of Aip, and D(At) of At. Plots of D(Aip) 
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FIG. 3: Lyapunov exponents vs. shell size e s for applied strain 
rate e a = 5.56 x 10 -6 s -1 ; embedding dimension d — 5. 



are shown in Fig. 4 for appropriately chosen values of 
e a . Peaked distributions are seen in the chaotic regime 
(5.6 x 10 _6 s _1 to 1.4 x I0~ 4 s~ 1 ) indicating the existence 
of characteristic values. Now we turn to higher strain 
rates. On the basis of single crystal results, we expect 
to find SOC state. In SOC type dynamics, power law 
statistics arises when spatially extended driven systems 
naturally evolve to a critical state. On increasing i a , the 
distributions become broader and asymmetrical in the 
mid-region of Aip eventually leading to power law distri- 
butions (Fig. 4c). For the data set at e a = 5.6xl0 _3 s _1 , 
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Distributions \D(Aijj)\ for the data sets at (a) 
\ (b) 1.4 x lO^s -1 and (c) 5.56 x 10" 3 s 



the distribution has the form D(Aip) ~ Aip~ a over one 
order of magnitude in Aip with a ~ 1.5. Similarly, we 
find a power law distribution D(At) ~ At~ b for the du- 
ration of the bursts with b ~ 3.2 and for the conditional 
average < Aip > c ~ At x with x ~ 4.2. Thus, the scaling 
relation x(a — 1) + 1 = b that characterizes SOC dynam- 
ics is well satisfied. Similar results are also found at the 
highest strain rate e a = 1.4 x 10~ 2 s _1 . Thus, the region 
of SOC-type dynamics extends from e a — 5.6 x 10~ 3 s~ 1 
onward, coinciding with the region of type A bands. 

Recall that a multifractal analysis was used to quantify 
the broad distribution of length scales occurring in the 
transition region between localized and delocalized states 
in the Anderson transition |d|. In our case also, there 
is a broad distribution of time scales in the region sep- 
arating chaos and SOC (Fig. 4b). Thus, we anticipate 
that multifractal analysis will quantify this heterogene- 
ity [pi. Let N = N(5t) be the number of time inter- 
vals at with m points, required to cover [0, if]. Then, 
the normalized amplitude of the bursts in the ith inter- 
val St, is Pi(St) = Y0k=i ^im+h/ Sj=i 4>j, usually called 
the probability measure. A conventional fractal is ade- 



quately described by a single scaling exponent, Df, its 
fractal dimension, through N(5t) ~ 5t~ D f . However, 
heterogeneous sets require a continuum of scaling indices 
introduced via the probability Pi(5t) ~ St a , where a is 
the strength of the local singularity. Then, N(St) gen- 
eralizes to N a (St) ~ 8t~f( a \ where f(a) is the fractal 
dimension of the subset of intervals characterized by the 
exponent a fig] . The non-uniformity of the measure is 
captured by the range of multifractality 9 — a max — a m i n , 
where a m in and a max are the extreme values of a. Using 
another measure Hi(6t,q) = pl/^2jPj, where q is a real 
number, a and f(a) can be directly calculated through 

m(St,q)lnpi(5t) 
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E t Vi(8t,q)lnp,i(6t,q) 
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This canonical method was shown to be suitable for the 
analysis of short experimental data sets ]T^| . The spec- 
trum (a, /(a)) is shown in Fig. 5 for e a = 1.4 x 10 _3 .s _1 . 
The dependence of the multifractal range 9 on e a is shown 
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in Fig. 6. Also displayed are the different dynamical 
regimes together with the band types. 9 is seen to have 
relatively low values at both low and high strain rates. 
The small value of 9 in the chaotic regime is due to the 
sharp peak in D (Aip) (Fig. 4a), while that in the SOC 
region is due to the scaling nature of the distribution 
(Fig. 4c) (ll|. In contrast, a sharp peak observed at in- 
termediate i a clearly signals a transition in the nature of 
dynamics from chaotic B-type to SOC A-type bands. 
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FIG. 6: The multifractal range 8 vs. applied strain rate i a - 
Regions of chaotic type B and SOC type A bands are marked. 

Distinct dynamical regimes and spatio-temporal pat- 
terns are displayed in the above analysis. Chaos is asso- 
ciated with type B hopping bands at low applied strain 
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rates and SOC dynamics is present at high e a in the do- 
main of type A propagating bands. The crossover from 
chaotic to SOC dynamics is clearly signaled by a burst 
in multifractality. Such a diversity in dynamics poses a 
challenge for modeling the PLC effect. What needs to 
be understood is why the spatial correlations between 
bands and the complexity of the dynamics increase with 
increasing e a . The following arguments are quite gen- 
eral and encompass various mechanisms that have been 
proposed to explain the spatial coupling between bands 
(see @, [§)• A band of localized deformation induces 
long range stresses that favor strain uniformization in 
its neighborhood. However, during jerky flow (i.e., even 
in the presence of solute atoms) relaxation of these in- 
ternal stresses is observed It occurs by thermally 
activated rearrangements of dislocation microstructure. 
Thus, both the magnitude and spatial extent of internal 
stresses are expected to decrease with time once a band is 
formed. Let t r denote this characteristic relaxation time. 

At very low strain rates, the reloading time between 
two successive drops, t l: is very large and ti >> t r . The 
internal stresses are fully relaxed and no spatial corre- 
lation between bands is observed. This corresponds to 
the domain of type C bands, which are not studied here. 
The absence of spatial correlation should make it easier 
to build a model based on the temporal mechanism of the 
instability, namely the repeated pinning and unpinning 
of dislocations in the field of solute atoms. 

At larger strain rates, the reloading time ti decreases 
and becomes commensurate with the relaxation time t r . 
Internal stresses are not totally relaxed and favor the 
formation of new bands nearby the previous ones, hence 
the hopping character of the associated type B bands. 



At high e a , ti « t r very little plastic relaxation occurs 
between the stress drops. Further, the reloading rates be- 
come increasingly a significant fraction of the unloading 
rates. Thus, the stresses felt by dislocations always re- 
main close to the threshold for unpinning from the solute 
atoms. New bands can be formed in the field of unrelaxed 
internal stresses and recurrent partial events can overlap 
resulting in a hierarchy of length scales. This leads to 
both SOC dynamics and type A propagating bands. 

Models incorporating the basic instability mechanism 
and a spatial coupling [Q, ^ seem to reproduce the var- 
ious types of bands and the corresponding stress drop 
distributions. This suggests that they have the right ba- 
sis for reproducing the corresponding dynamical regimes, 
chaos and SOC, as well as the crossover between them. 

In summary, a connection is made between the spatial 
patterns and the dynamical regimes associated with dif- 
ferent strain rates in the jerky flow of polycrystals. A 
qualitative interpretation is proposed in terms of com- 
peting mechanisms that operate at different scales, local 
and global, and involve non-local effects. The crossover 
between the hopping bands associated with chaos and the 
propagating bands identified with the SOC state is well 
detected by a surge in the multifractal behavior of the dis- 
tribution of plastic events. To the best of our knowledge, 
this is first example where a crossover from localized to 
propagating states has been detected using a multifractal 
analysis based on purely experimental signals. 
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